Krylov iterative methods and synthetic acceleration for transport in binary statistical media
نویسندگان
چکیده
Iterative solution of the Levermore-Pomraning equations for transport in binary statistical mixtures can be extremely slow in certain limits. We propose an iterative method that improves convergence by utilizing a combination of inner iterations, synthetic acceleration schemes and Krylov iterative methods. Spectral analysis and numerical results show that our new scheme outperforms simpler iterative methods for problems in either the diffusion or the atomic mix regimes, or both.
منابع مشابه
New variants of the global Krylov type methods for linear systems with multiple right-hand sides arising in elliptic PDEs
In this paper, we present new variants of global bi-conjugate gradient (Gl-BiCG) and global bi-conjugate residual (Gl-BiCR) methods for solving nonsymmetric linear systems with multiple right-hand sides. These methods are based on global oblique projections of the initial residual onto a matrix Krylov subspace. It is shown that these new algorithms converge faster and more smoothly than the Gl-...
متن کاملPreconditioned Krylov Subspace Methods for Transport Equations
Transport equations have many important applications. Because the equations are based on highly non-normal operators, they present diiculties in numerical computations. The iterative methods have been shown to be one of eecient numerical methods to solve transport equations. However, because of the nature of transport problems, convergence of iterative methods tends to slow for many important p...
متن کاملAcceleration of GPU-based Krylov solvers via data transfer reduction
Krylov subspace iterative solvers are often the method of choice when solving large sparse linear systems. At the same time, hardware accelerators such as graphics processing units continue to offer significant floating point performance gains for matrix and vector computations through easy-to-use libraries of computational kernels. However, as these libraries are usually composed of a well opt...
متن کاملQuasi-Newton acceleration of ILU preconditioners for nonlinear two-phase flow equations in porous media
In this work, preconditioners for the iterative solution by Krylov methods of the linear systems arising at each Newton iteration are studied. The preconditioner is defined by means of a Broyden-type rank-one update of a given initial preconditioner, at each nonlinear iteration, as described in [5] where convergence properties of the scheme are theoretically proved. This acceleration is employe...
متن کاملA Semi-Lagrangian Two-Level Preconditioned Newton-Krylov Solver for Constrained Diffeomorphic Image Registration
We propose an efficient numerical algorithm for the solution of diffeomorphic image registration problems. We use a variational formulation constrained by a partial differential equation (PDE), where the constraints are a scalar transport equation. We use a pseudospectral discretization in space and second-order accurate semi-Lagrangian time stepping scheme for the transport equations. We solve...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- J. Comput. Physics
دوره 228 شماره
صفحات -
تاریخ انتشار 2009